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Abstract 

We consider the steady Swift - Hohenberg partial differential equation. It is a one- 
parameter family of PDE on the plane, modeling for example Rayleigh - Benard con- 
vection. For values of the parameter near its critical value, we look for small solutions, 
quasiperiodic in all directions of the plane and which are invariant under rotations of 
angle 7r/q, q > 4. We solve an unusual small divisor problem, and prove the existence 
of solutions for small parameter values. 

Contents 

1 Introduction 2 

1.1 Main result and sketch of the proof [3 

2 Analysis of the main part of the differential 5 

2.1 Setting 5 

2.2 Formal computation 7 

2.3 Formulation of the problem 8 

2.4 Splitting of the space and first reduction of the problem [S 

2.5 Reduction to the subspace E2 1C 

2.6 Structure of the reduced operator A e 16 

2.6.1 General structure - Invariant subspaces 16 

2.6.2 Computation of coefficients 7^ in Ai 21 

2.7 Eigenvalues of A e 22 

2.8 Inverse of C e in U s m 



3 Existence of the solution 24 



'University of Groningen, Department of Mathematics P.O. Box 800, 9700 AV Groningen, The Nether- 
lands. E-mail : braaksma@rug.nl 

^Institut Universitaire de France-Laboratoire J. -A. Dieudonne U.M.R. 7351, Universite de Nice - Sophia 
Antipolis, Pare Valrose 06108 Nice Cedex 02, France. E-mail : gerard.iooss@unice.fr 

"^CNRS-Laboratoire J. -A. Dieudonne U.M.R. 7351, Universite de Nice - Sophia Antipolis, Pare Valrose 
06108 Nice Cedex 02, France. E-mail : stolo@unice.fr. Ce travail a beneficie d'une aide de l'Agence 
Nationale de la Recherche portant la reference "ANR-10-BLAN 0102" 



1 



A Proof of Lemma [6] 
B Proof of Lemma 1161 



25 







1 Introduction 

In the present paper we study the existence of a special kind of stationnary solutions 
(i.e. independent of t), bifurcating from (i.e. tending towards zero when the parameter A 
tends towards 0), called quasipatterns of the 2-dimensional Swift-Hohenberg PDE 



and A is a parameter. These are two-dimensional patterns that have no translation sym- 
metry and are quasiperiodic in any spatial direction. 

Mathematical existence of quasipatterns is one of the outstanding problems in pattern 
formation theory. To our knowledge, hereafter is the first proof of existence of such quasipat- 
terns of a PDE. Quasipatterns were discovered in nonlinear pattern- forming systems in the 
Faraday wave experiment |BCM92 l [EF94J, in which a layer of fluid is subjected to vertical 
oscillations. Since their discovery, they have also been found in nonlinear optical systems, 
shaken convection and in liquid crystals (see references in |AG12| ) . In spite of the lack 
of translation symmetry (in contrast to periodic patterns), the solutions are 7r/g-rotation 
invariant for some integer q (most often observed, 2q is 8, 10 or 12). 

In many of these experiments, the domain is large compared to the size of the pattern, 
and the boundaries appear to have little effect. Furthermore, the pattern is usually formed 
in two directions {x\ and X2), while the third direction (z) plays little role. Mathematical 
models of the experiments are therefore often posed with two unbounded directions, and 
the basic symmetry of the problem is E(2), the Euclidean group of rotations, translations 
and reflections of the (x\,X2) plane. 

The above model equation is the simplest pattern- forming PDE, and is extremely suc- 
cessful for describing primary bifurcations (the first symmetry breaking) of hydrodynamical 
instability problems such as the Rayleigh - Benard convection. Its essential properties are 
that 

i) the system is invariant under the group E(2); 

ii) the instability occurs for a certain critical value of the parameter (here A = 0) for 
which critical modes are given by wave vectors sitting on a circle of non zero radius (here 
the unit circle); 

iii) the linear part is selfadjoint and contains the main derivatives. 

In contrast to periodic patterns, quasipatterns do not fit into any spatially periodic 
domain and have Fourier expansions with wavevectors that live on a quasilattice (defined 
below). At the onset of pattern formation, the critical modes have zero growth rate but 
there are other modes on the quasilattice that have growth rates arbitrarily close to zero, 
and techniques that are used for periodic patterns cannot be applied. These small growth 
rates appear as small divisors, as seen below, and correspond at criticality (A = 0) to the 





where u is the unknown real- valued function on some subset of IR + xK , A: 




) 
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Figure 1: Example 8- fold quasipattern after |IR10] . This is an approximate solution of the 
steady Swift-Hohenberg equation ([2]) with A = 0.1, computed by using Newton iteration 
to find an equilibrium solution truncated to wavenumbers satisfying |k| < y/b and to the 
quasilattice T27 obtained with iVk < 27. 

fact that for the linearized operator at the origin (here —(1 + A) 2 ), the eigenvalue is 
not isolated in the spectrum, appearing as part of the continuous spectrum. 

If a formal computation in powers of vA is performed in this case without regard to its 
validity, this results in a possibly divergent power series in the parameter, and this approach 
does not lead to the existence of quasipattern solutions, but instead to approximate solutions 
up to an exponentially small error fIRlOj. 

In this work, we prove the existence of quasipattern solutions of the steady 
Swift-Hohenberg equations. Our result rests on the article [IR10| by G. Iooss and A.M. 
Rucklidge which settle the mathematical foundation of the problem such as the formulation 
of suitable functions spaces. We refer to the articles of Rucklidge [RR03 ; RS07] and Iooss- 
Rucklidge[IR10] for physical motivation as well as for the bibliography. 

1.1 Main result and sketch of the proof 

The problem is to find a special kind of solutions defined on M? of the steady Swift- 
Hohenberg equation 

(1 + A) 2 U - \U + U 3 = 0. (2) 

The parameter A is supposed to be real and small in absolute value. The solutions we are 
interested in should tend towards zero as the parameter goes to zero. 

We study equation ([2]) for A > 0. Namely, let Q = 2q be an even integer and let 
kj = exp , j = 1, . . . , 2q be the 2q unit vectors of the plane, identified with roots 

of unity. Let T be the set of linear combinations of vectors kj with nonnegative integer 
coefficients. We look for the existence of a (nonzero) 7r/(/-rotation invariant solution of the 
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form 

U(x) = Y^u {ic) e ik - x 
ker 

which belongs to a "Sobolev" like space 7i s , s > : 

\\Uf s :=^|n( k )| 2 (l + iV2r <+oo. 
ker 

The natural number denotes the minimal length of the linear combinations of the kj 's 
needed to reach k. 

We then show that such a solution exists indeed, for small positive parameters A: 

Theorem 1. For any q > 4 and any s > q/2, there exists Xq > 0, such that the steady 
Swift- Hohenberg equation for < A < Aq, admits a quasipattern solution in % s , invariant 
under rotations of angle ir/q. Moreover the asymptotic expansion of this solution is given 
by the formal series given in flRlO]. 

One of the main difficulty is that the linearized operator at W = 0, has an unbounded 
inverse. Indeed, it is easy to show that the eigenvalues of (1 + A) 2 in T~L S are (1 — |k| 2 ) 2 
where k S T. These numbers accumulate in 0. It creates a small divisor problem, such 
that if A > nothing can be said a priori about (1 + A) 2 — Aid. We use the first terms 
of the asymptotic expansion of the solution and change the unknown as U = U € + e 4 W 
and A = e 2 (A2 + A4e 2 ) for some well chosen U e and positive A2. Let C e be the linear part 
at W = of the nonlinear equation so obtained. For e = 0, the operator Cq = (1 + A) 2 
is a positive selfadjoint operator in Ti s . It is bounded from T-L s +i into rl s , but it is not 
Fredholm, since its range is not closed. Its spectrum is an essential spectrum filling the half 
line [0,oo). The set of eigenvalues is dense in the spectrum. The linear operator C t is the 
sum of Co and a bounded operator (multiplication by a small function 0(e 2 )) selfadjoint in 
%q. If the range of C e were closed, a usual way to estimate the inverse of the selfadjoint 
operator C t in Hq, would be to estimate the distance from to its numerical range (see 
|Kat95| ) (containing the spectrum). Such estimate as 

(C t U,U) >ce 2 \\U\\ 2 

for a certain constant c > 0, cannot be proved here. So we need to study the linear operator 
in more details. 

We show that there exists an orthogonal decomposition (depending on e) of the space 
Us = Eq © Ei © E2, s > q/2, such that the solution of the equation CJJ = f in rl s can 

be computed and estimated from its ^-component C/ 2 . The latter is solution of a linear 

(2) 

equation C e U2 = f ■ The main part (with respect to powers of e) of that operator is 
an operator A e . When restricted to rotation invariant elements, it has a block-diagonal 
structure of fixed finite dimensional blocks. Then, it is possible to estimate all eigenvalues 
of these selfadjoint blocks. These eigenvalues have the form 

(|k| 2 -l) 2 + 3e 2 + 0(e 4 ) 

for keT. One of the main feature is that they do not accumulate at the origin, and 

despite of the small divisor problem arising for e = 0, we are able to give an upper bound 
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in Ho of the inverse of A e of order 1/A for nonzero A. Then we extend this estimate in H. s 
for s > q/2. Finally, we use a variant of the implicit function theorem to conclude to the 
existence of quasipatterns solutions. 

Remark 2. If the coefficient (3) of e 2 were negative, then the operator could have again 
small divisors. Then the proof of the existence would have been much more involved. At 
that point, solving the nonlinear problem in U2 would have required the use of a version of 
Nash-Moser theorem such as those developed by J. Bourgain, W. Craig, M. Berti and al. 
(see for instance ]Bou9b\ I CraO(A \BBP1C\ BerOT^ ). Their main feature is the use of the good 



separation property of the "singular sites" of the main linear operator. Indeed, we can show 
that the operator (1 + A) 2 does have this property. 

Remark 3. For A < the solution U = is isolated in an open ball of radius as it 

is easily deduced from the following estimate 

HKl + Aja-Art^lAl- 1 

which holds in He. 



2 Analysis of the main part of the differential 
2.1 Setting 

In this section, we recall and improve some of the properties of the function spaces we use, 
as defined in |IR10j . 

Let Q = 2q, q > 4 being an integer. Let us define the unit wave vectors (identifying C 
with M 2 ) 

■ j-i 

kj := e n 1 , j = 1,. . . ,2q. (3) 

We define the quasilattice r C I 2 to be the set of points spanned by (nonnegative) 
integer linear combinations of the kj 's : 

2q 

k m = ^m j k i , m = (mi,...,m 2 ,)eN 2? . (4) 
3=1 

We have kj = —kj +q . Hence, we can write 

9 

km = ^ ^ m j^j 1 

i=i 

where, m'j := rrij — rrij +q belongs to Z. Thus, 

<? 

|k m | 2 = m^m'j < ki,kj > . 

We then define, for any m E N 2q and k G T, 

|m| := rrij, := min{|m| : k = k m }. 

j 

We have 
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Lemma 4. [IRi($ [Lemma 4-1] For any kef, we have : 



iVk+k' < iVk + N v , N_ k = N k 



(5) 



|k| < iV k 

card{k | N k = N} < ^{q)^' 1 
for some constant ci(q) depending only on q. 
As in [IR10] . we use function spaces defined as 



J W = Y J W^e ikx ; \\W\\ 2 S = ^(1 + N^YlW^l 2 < oo 
I ker ker 



which are Hilbert spaces with the scalar product 

(w,v)s = ^(i + iv k 2 )^ (k) F 

ker 



(k) 



Lemma 5. For s > q/2, for any U E % s and any V G T~Lq, we have 

\\UV\\o < c s \\U\\s\\V\\ 

for a certain constant c s > 0. 

Proof. Using Cauchy-Schwarz inequality, we have 



ic^iig < E 



ker 



V [/ (m V (k_m) 



mer 



(6) 



(7) 



(8) 



0) 



< E(Ei c/(m) i 2 ( 1 + iV -) s ) f E w^'W+N^r 8 ) 

ker Vmer / Vm'er / 

< ii^ii^Ei y(k " m,) i 2 )^ 

Vker / 

where R s := £ m , gr (l + ^m')^- This last sum converges if s > q/2. Indeed, according to 
[IRT0] [(24)]. card{k G T \ N k = N} < c(q)N q ~ l for some constant c{q). Hence \\UV\\l < 

\qFLs- 



\\U\\ 2 s \\V\\lR s . | 

Lemma 6. (Moser-Nirenberg type inequality) Let s, s' > q/2 and let U, V £ l~L s r)H s '. Then, 

\\UV\\ S < C(s,s')(\\U\\ s \\V\\ s , + ||J7|| S '||F|| S ) (10) 

for some positive constant C(s, s') that depends only on s and s' . For i > and s > l + q/2, 
% s is continuously embedded into C 

We postpone the proof to the appendix. 
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2.2 Formal computation 

Let us look for formal solutions of the steady Swift-Hohenberg equation 

XU - (1 + A) 2 U - U 3 = 0, (11) 

We characterise the functions of interest by their Fourier coefficients on the quasilattice T 
generated by the 2q equally spaced unit vectors k, (see Q): 

[/(x) = ^ U ( k )e ik ' x , x=(i 1 ,s 2 )el 2 . 
ker 

We seek a non trivial solution, bifurcating from 0, parameterized by e, and which is 
invariant under rotations by n/q. As it is shown for example in |IR10j . a formal computation 
with identification of orders in e leads to 

U(xi,x 2 ) = eu (xi,x 2 ) + e 3 ui(xi,x 2 ) + ■ ■ ■ A = e 2 A 2 + e 4 A 4 + . . . (12) 

and gives at order 0(e) 

= (l + A) 2 -u - (13) 
We take as our basic solution a quasipattern that is invariant under rotations by n/q: 

2q 

uo = X>* J ' X - (14) 

At order 0(e 3 ) we have 

A 2 u -ujj = (1 + A)V. (15) 

In order to solve this equation for u%, we must impose a solvability condition, namely that 
the coefficients of e* kj " x , j = 1, ...,2q on the left hand side of this equation must be zero. 
Because of the invariance under rotations by ir/q, it is sufficient to cancel the coefficient of 
giki-x p or ^ e com putation of the coefficient, we need the following property 
Property: If we have 

kj + k ; + k, r + k s = for j, l,r,s G {1, 2q} 

then either kj + kj = 0, or kj + k r , or \tj + k s = (there are two pairs of opposite unit 
vectors). 

Proof. Since there are 4 unit vectors on the unit circle, we can assume without restriction, 
that kj and k; make an angle 26 < tt/2. Then |kj + kj| = 2cos# > y2. It results that 
|k r + k s | = 2cos^ with k r and k s symmetric with respect to the direction of the bissectrix 
of (kj,kj), making the same angle as kj and k; with the bissectrix. So {k r , k s } is the 
symmetric with respect to of {kj, k^}. | 

This yields 

A 2 = 3(2g-1) (16) 
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which is strictly positive. Moreover 

u x = ° kelk ' x > ( 17 ) 

ker,iv k ^i,7V k <3 

3 

"3k, = -1/64, o^kj+k, = ~ (1 _ | 2k +kf | 2 )2 ' k i + k ^0, 

«k, +k , + k,. = - (1 _| k . + 6 k; + kr |2 )2 >i^^^^ 

+ k, ^ 0, kj + k r / 0, k r + k, / 0. 

We notice that for any k, ak < in u%. At order C(e 5 ) we have 

A4110 + A 2 ni - 3tigni = (1 + A) 2 u 2 . (18) 

The solvability condition gives A4 equal to the coefficient of e* kl ' x in 3uqU\. From the 
expression (fTTl) of ui 

A 4 = a k 

k J +k i +k=ki, AT k =3 

where all coefficients ak are negative, it results that 

A 4 < 0. (19) 

2.3 Formulation of the problem 

Let us define the new unknown function W in rewriting (|12p as: 

U = U e + e 4 W 

U e = eu + eV + e 5 u 2 (20) 
A e = e 2 A 2 + e 4 A 4 

where no, ui, u 2 , A 2 , A4, are as above. Given a particular (small) positive value of A, we 
get e 2 by the implicit function theorem, and since A 2 > 0, we obtain a unique positive e. 
All the corrections are in W. The aim is to show that the quasi-periodic function W exists 
and is small as e tends towards 0. By construction we have 

\JJ e - (1 + A) 2 U e - U 3 =: -e 7 f e 

where f t is quasi-periodic, of order 0(1) with a finite expansion, and is function of e 2 . After 
substituting (j20|) into the PDE (fTT]) . we obtain an equation of the form 

F(e,W) = 0, 

with 

F(e, W) =: C e W + e 3 f e + 3e 4 U e W 2 + e 8 W 3 , (21) 

where 

C € = (l + A) 2 -A e + 3C/ e 2 =:£ e + e 6 P £ , (22) 
£ e = (l + A) 2 + e 2 a + e 4 6, (23) 

a = 3uq — A 2 , b = GuqUi — A4, 
V e = 6u u 2 + 3(ui + e 2 « 2 ) 2 . 
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Remark 7. The degree of truncation, that is the degree in e in the expansion of U e , is 
chosen so that the power of e in front of both W 2 and f e are greater than 2. This is crucial 
for the very last step of the proof. 

It is clear that the operator V e is an operator bounded in any T-L r , r > uniformly 
bounded in e, for e < eo- 

A nice property of the operator £ e is that the averages ao and bo of a and b are strictly 
> 0. Indeed, we have for any q 

ao = 3, bo = — A4 > 0. 

For ao this results from (I16p and a simple examination of Uq the average of which is 2q, and 
for 60 we observe that the average of uou± is 0, due to the form of u\. 

Assume that we could prove that £~ is 0(e~ 2 ). Then, provided that we are in H s , s > 
q/2 which is a Banach algebra, we should get from (|2ip 

W = 0(e) + 0(e 2 \\W\\ 2 ) 

The standard implicit function theorem then would allow to conclude and to get W = 0(e). 

In fact, it is not expected that the operator C e has a bounded inverse, due to the small 
divisor problem mentioned in section 1.1. Notice in particular that it is shown in [IRIOJ 
that there exists c > such that for any k G T\{kj; j = 1,2, ...,2g} 

llkl 2 -!!^^, (24) 

where Iq + 1 is the order of the algebraic number uj = 2cos7r/g. This estimate is similar to 
the Siegel's diophantine condition for linearization of vector fields |Arn80| . 

This lower bound shows that the inverse of £0 on the orthogonal complement of its 
kernel is an unbounded operator in H s , only bounded from ~H S to H s -4i - In other words, 
belongs to the continuous spectrum of £0 an d the main difficulty to be solved below is to 
find a bound for the inverse C~ l for small values of e. Notice that C e is selfadjoint in TLo 
but not in T~L S for s > 0. It is tempting to work on its small (real) eigenvalues to obtain 
a bound of its inverse. However, we are in infinite dimensions, so the spectrum does not 
contain only eigenvalues, and an option would be to truncate the space to functions with 
finite Fourier expansions (with k such that < N) . Since our method consists in reducing 
the study to an operator in a smaller space, it is preferable to use the eigenvalues later, on 
the reduced operator. 

2.4 Splitting of the space and first reduction of the problem 

Let us split the space H s into three mutually othogonal (in any % s ) subspaces. We define 



Eo 



Ei 



Eo 



J W = ^ (k) e ikx G U s \ ||k| 2 - 1| 2 > 5 2 , and |k - k,-| > 6 1} ,j 6 {1, 2q} 1 
I ker J 

J w = ^ (k) e ik x G U 3 \ k G <ti 1 , 
I ker J 

J W = Y VF (k) e ikx G n s ;3j G {1, ..2q} such that k G a 2 j \ , 
I ker J 
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Figure 2: Division of the Fourier spectrum into o"o, <J\ and a 2 . 

where (see figure [5]) 

a x = {k G r; ||k| 2 - 1| 2 < 5 2 and |k-k,| > Si = V35,j G {l,2g}}, 

o"2,i = {k G T; |k - kf| < <ji}, a 2 = U 2q =l a 2 j. 

In figure [2] the annulus centered at the unit circle has a thickness 25 and the little discs 
should have a radius <5i = V35, such that the intersection of the shaded area with the shifted 
one, centered at the point (2ki,0) is, for 5 small enough, reduced to the disc centered at kx- 
In the sequel, we choose 5 = Ce 1 / 2 with C large enough and 

hence 25 + 5 2 < 5\ is verified for e 1//2 < 1/C and the intersection o~\ D {o\ + 2ki} is empty 
(hint: solve 5\ = (1 + 5) 2 — 1 for intersecting the circle centered in 0, of radius 1 + 5 with 
the line of abscissa 1 parallel to y axis). 
This leads to (see figure [3|) 

o-i n {(7i + kj + k z , j,l = 1, ..2q, kj + k ; / 0} = 0. 

The subspaces Ei are closed as intersections of closed subspaces (kernel of certain coef- 
ficients, still continuous functionnals here) and we have the orthogonal decomposition 

U s = E ®E 1 ®E 2 . 

The orthogonal projections associated with this decomposition are denoted by Pq, P\, P 2 . We 
also notice that the multiplication operator by a function having a finite Fourier expansion 
with wave vectors in T is a bounded linear operator in H r for any r > 0. Indeed a finite 
Fourier expansion belongs to % s for any s, and Lemmas [5] and [6] apply. 

2.5 Reduction to the subspace E 2 

The aim here is to solve with respect to U G % s the equation 

£ e U = f, (25) 
where C € is defined in (|22j) and where / G 7i s is given. 
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Figure 3: Empty intersections of u\ with o\ shifted by kj + 
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Remark 8. We denote by Co the operator £0,0 = £-o,w = (1 + A) 2 = £o- For being more 
explicit in formulae, we denote by PjC^Pj the restriction of Cq to the invariant subspace 
Ej,j = 0,1,2. 

We decompose this equation on the subspaces Eq, E\, E 2 , which gives, after noticing 
that the subspaces Eq, E\, E 2 are invariant under C , and that £o|.Eu A) I £2 are bounded 
operators, and PjC^Po = for j € {1, 2}, 

P C e U + P {(e 2 a + e 4 b + e 6 V e )(U 1 + U 2 )} = f (26) 
P 1 C e U l +P l {(e 2 a + e 4 b + e 6 V e )(U + U 2 )} = h (27) 
P 2 C e U 2 + P 2 {(e 2 a + e 4 b + e 6 V e )(U + U l )} = f 2 (28) 

where f j = P j f,j = 0,1, 2. 
We have the following 

Lemma 9. Let fix S > and choose s such that < s < S. Then, there exists €q > 0, such 
that choosing C large enough in the definition of 5, e < eo, the equation CJJ = f in H s 
reduces to 

Cf)u 2 = Q(e)f 

where U 2 = P 2 U, 

4 2) = A e + e 4 R e 
A t U 2 = P 2 (C + e 2 a)U 2 

K e U 2 = P 2 {bU 2 -a{P Q C Q P Q y l Pv{aU 2 )} + 0{e 2 )U 2 , (29) 

Here 0(e 2 ) denotes a bounded linear operator in T{ s which norm is bounded by ce 2 . The 
components Uq := PqU and U\ := P±U are functions of U 2 and f and satisfy the following 
inequalities : 

\\Uo\\s < ce 2 ||C/ 2 || s + ^||(Po+Pi)/||, 
\WA\s < ce 4 \\U 2 \\ s + ^(ePo + P^fWs 

where a = a — oq, for a certain c > 0, only depending on s, and R e and Q(e) are bounded 
linear operators in H s , depending smoothly on e. 

We need the following Lemma 

Lemma 10. For e small enough, we have 

aUx^Eo, Pi{aU 2 ) = 0, Pi(W7 2 ) = 0, (30) 

where a is the oscillating part of a : 

a = a — oq. 

Moreover, the Fourier spectra of Po(aU 2 ) and Po(bU 2 ) are at a distance of order 1 of the 
unit circle. 
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By distance of order 1 of the circle, we mean a strictly positive distance as e tends to 0, 
in the decomposition into subspaces Ej. 

Proof. We need to prove that i) for k £ o~\, then k + k m ^ a\ U o 2 for k m 7^ 0, |m| = 2, 
and ii) for k £ o~2, then k + k m £ cti with |m| = 2 or 4. 

For showing this, let us first observe that the intersections of the unit circle with all 
circles of radius 1, centered at k m , with k m ^ 0, |m| = 2, are exactly the 2q points k r , 
r = 1, 2q. Let us define the region do (shaded region in figure [2] (a)) defined by 

a = cri Uj- e { lr . i2g } a 2 j 

which is the union of the annulus o\ and the discs Uj 6 {i 2<j} cr 2j - Now consider the inter- 
section of do with the union of shifted analogue annuli Ok m 

a k m = {k + k m ; k £ o , k m ^ 0, |m| = 2} 

centered at all k m ,such that |m| = 2 (see figure [3|). It is clear that for any given q, and e 
small enough, the little discs of radius 5i are such that the intersection 00 H {U| m | =2 cik m } is 
exactly the union of the little closed discs centered at each k,. It results that for U £ E%®E 2 , 
the product all for which the corresponding wave vectors belong to some dk m , has a zero 
projection on E\. This proves that 

P 1 (aU 1 ) = 0, Pi(aU 2 ) = 0, 

which implies 

P l (aU 2 ) = 0. 

It is also clear that 

a n {k + k m ; k £ cti, k m / 0, |m| = 2} = 0, 
which moreover implies that 

P 2 {aU 1 ) = 0, 

and (|30p is proved for the part concerning a. Now observe that we have |k,- + k m | 7^ 1 except 
when kj + k m = k r for some r £ (1, 2q). Since we only consider the finite number of cases 
|m| = 2 or 4, it is clear that in choosing e small enough (i.e. 5± small enough), then for 
k £ a 2 , k + k m ^ o\. It results in particular that 

P 1 (bU 2 ) = 0. 

The last assertion of Lemma [10] results from the fact that |k + k m | 7^ 1 for the Fourier 
spectrum of terms £ ^0(0^2) and ^0(^2), with a distance to the unit circle equivalent to 
||kj + k m | — 1| when it is not zero. | 

Proof, of Lemma [9] : We know by construction that 

\\(P C P )-%<-±-, 
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Hence, we have 

(PoCePoy 1 = [1 + {PoLoPotH^a + e% + e & V e )]- l {P Q C Q P Q y\ (31) 
and the estimate 

\\e 2 a + e 4 b + e 6 V e \\ s < c(s)e 2 
leads for e small enough (s < S) to 

(PoC.Pq)- 1 = [I-(P C P )- 1 (e 2 a + e 4 b + e (i r € ) + 



+{(PoC Por\e 2 a + e A b + e 6 V e )} 2 + 0( e 3 )](P £ Po) -1 



with a convergent Neumann power series in the bracket, as soon as ^f)r < 1, which holds 



for e small enough (s < S). The first consequence is 

1 



(P C e P )-% < 



e(C 2 -c(s)e)' 
Notice that 

Uq = -(P^P^P^a + e% + e 6 P e )(£/i + U 2 )} + (P £ e P ) -1 /o (32) 

The last property of Lemma [10] and (f3Tj) imply that in (f32|) we have, for e G (0, eo) and for 

any s 

\\(P C Po)- 1 Po(aU 2 )\\s < d(s)\\U 2 \\ s , 
WiPoCoPor'PoibUz^s < d(s)\\U 2 \\ s , 

\\(P C Po)- 1 a(P CoP )- 1 Po(aU 2 )\\ s < d(s)\\U 2 \\ s . 

It results that 

-(Po^Po)" 1 Po{(e 2 « + e 4 6 + e 6 P e )f/ 2 } = -^{P^P^P^aU^ + 

+0(e 4 \\U 2 \\ s ), 

hence 

U Q = Q Q ,i{e)U x + Qo, 2 (e)C/ 2 + (iW%) _1 /o (33) 

with 

Qoj(e) =: -(PoAPo)- 1 Po(e 2 a + e 4 6 + e 6 P e )P i; j = 1, 2, 
and for s < S, e £ (0,€q(S)) the following estimates hold: 

Qo, 2 (e)C/ 2 = -e 2 (P £ Po)" 1 Po(a^ 2 ) + e 4 Qo. 2 (e)^ 2 , 



\\Qo,iWl\\s < ^Z^ s)£ \\Ui\\s < ^eims, (34) 
||^, 2 (e)f/ 2 || s < c^WUzWs, 

(PoCePo)^ 1 fo\\ s < H/olU- 
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It results from Lemma [TO] that equation ([27|) leads to 

Pl(£o + e 2 a + e 4 6 + e % V,)U x + Pi{(e 2 a + e 4 b + e 6 V 6 )U } + Pie 6 V e U 2 = fx (35) 

and since ao = 3 > and the operator P\CqP\ is positive, we can invert the operator 
-Pi (A) + e 2 ao)Pi with the estimate 

\\{P 1 (c + e 2 a )p 1 }- 1 \\<^. 

Now replacing Uq by its expression (|33p into equation (j35h . we introduce an operator acting 
on U\ of the form 

Pi(e 4 5 + e 6 P e )Pi + Pi(e 2 « + e% + e 6 P e )Q ,i(e) 

which is bounded by 0(e 3 ), perturbing Pi(£o + e 2 oo)Pi the inverse of which is bounded by 
l/3e 2 . It results that, for e small enough, the operator acting on U\ has a bounded inverse, 
with 

\\{Pi[C e + (e 2 a + e 4 6 + e^QoA^PiVX < \- 

Moreover by Lemma [TOl we have Pi ^{PqCqPq)^ 1 Po{aU 2 )] = 0, hence it results that there 
are bounded linear operators 

:E 2 ^E , BP :E 2 ^E 1 , 

such that 

U x = e 4 BW(e)U 2 + -Q^(e)f + ^Q^ 1 \e)f 1 , (36) 

U = e 2 B^(e)U 2 + -Q^{e)h + -Q^(e)h, (37) 

e e 

with the estimates 

||^')(e)|| s < c 2 (s),j = 0,1, 
||Q ( ^(e)IU < c 2 (s), i,j = 0,l, 
uniform in e G (0, e ) and in bounded in % s . Moreover, as e — )• 

fj(0)( £ ) = 4 0) +O(6 2 ), 

^ 0) C/ 2 ~ -(Po£oPo) _1 Po(^2). 

Equation ([28]) now reads 

P 2 £ £ C/ 2 + P 2 {(e 2 a + e 4 6 + e 6 V e )U } + P 2 {(e 4 6 + e 6 7> £ )E/i} = / 2 

and replacing Uo and £7i by their expressions ([3"T]) . ([3T>]) in function of C/ 2 leads to 

P 2 £ e *7 2 + P 2 {(e 2 a + e 4 6 + e 6 V t )e 2 B^ (e)U 2 + 
+P 2 {(e 4 6 + e 6 P e )e%( 1 )(e)C/ 2 } 

= /2-P2((e 2 a + 6 4 6 + e 6 n)^Q (0 ' 0) (6)/o + ^S (0 ' 1) (e)/i]) + 

-P 2 ( (e 4 6 + e 6 n)[ I Q (i,o) (e)/o + ^ g (i,i) (e)/l] ^ ; 
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this gives (see the definition of A e in Lemma [9]) 

A £ U 2 + e 4 R e U 2 = Q(e)f 

with the announced properties for bounded operators R e and Q(e) in H s . Lemma [9] is 
proved. | 

2.6 Structure of the reduced operator A e 
2.6.1 General structure - Invariant subspaces 

We study in this section the structure of the operator A e defined by (|29p . We first observe 
that we deal with functions U which are invariant under the rotation Rn of the plane. So, 
let us define a new subspace of H s for such functions: 



J 2 

This implies immediately that 



E { 2 S) = {U e E 2 ; U(Rex) = f/(x)}. 



UM = U {R ^\ (38) 



(s) 

and for any U E E^ we have the following decomposition 

j = l,...,2q 

where 

U 2d (x)= £ C/ (k) e lk - x = C/ 2 ,i(i?(w).x). (39) 

(S) 

It results that any U £ E 2 may be written as 

C/(x) = ]T t/ 2 .l(fi(l-j)xX). 

Let us notice that in the little disc o"2,i we have 

£j 2 ,i 9 k = ki + k', |k'| < 5i = e 1 / 4 v / 3C, 

and let decompose the discs a 2) i into 2q equal sectors k^ + S TO; m = 1, 2q such that 

, / ^ /, c- ,/ r (m — l)-7r 7r (m — l)ir ir , ] , 

k € T; k < <5i, arg k G - — , ^ L + - . 40 

g 2q q 2q J 

For any k' 6 S;, we define the set of 2q spectral points ajj by 



a 



(0 



/ = {k = k j +k / Ga2, j ; j = l,...,2g}. 
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The subspace of % s associated with cr$ is denoted by E^ k , , so that any U G E^p may be 



written as 



= E E t/^ +k ')e l(kj+ ^ k0 ' X 

j'=l,...,2gki+k'6o- 2 ,i 

= E E E c/ (kl+i '^ k ' ) e i ^ +k ')- x , 

i=l,...,2 ? k'eS ( j=l,...,2q 



hence 



= E E c/a,k ' ) ' c/(i,k ' ) G 

j=l,...,2g j=l,-,2g 

and any f/ € is completely determined by the set of 2q— dimensional [/^ 1,k '- > G -E^k'' 
k' G Si, identified with the set of components 

{C /(k,+k) =[/ } J = l,...,2g}. 

Indeed we have 

^ k ')(x) = C/ (m ' R f k ' ) ( J R f x), 

hence 

^(x) = U {1 ' R ^'\ R{1 _ ir _x), k' G S, 

(/,R (i _ 1)7r k') (ki+fl ^j^ k') i(kj+fi (^i^ k , )-x 

?7 (x) = u ~^r~ e , k G Si, 

i=l,...,2g 

where 

^• k ') (x)= £ ^^^(k.+k')^ k , eSi) (41) 

j=l,...,2q 

and we observe that the coordinates of U » , k G Si, correspond to those shifted 
of L^ 1 '* 4 ). Moreover we have 

2q 

C/ ( X ) = E E ^ (1 ' k,) (%-0f x )- ( 42 ) 

Z=l k'GSi 

From now on, we denote by E 2 y the previously defined 2q-dimensional subspace E^,. 

Looking at the form of the operator A e we see that the wave vector k of U is shifted by 
k m , |m| = 2 at order e 2 , and |m| = 4 at order e 4 . Now, we observe that for a fixed finite 
|m|, if the combination 

k m - (ki - kj) 
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is not 0, then it has a minimal length of order 1 as e tends to 0. It results that for 
k = ki + k' G (72,i, and 1 = kj + 1' G a 2 j with 

k m - (k - 1) ^ 

then, for e small enough 

k m -(k-i) = e>(i). 

It results that the only possibility for going from 1 G (J2,j to <7 2; i is to add k m = ki +k J+(J = 
ki — kj. It results that the system 

k m + 1 = k, 1 = kj + 1', k = k, + k' 

has the only solution 

1' = k'. (43) 

It should be clear that k m comes from terms with many possible combinations, not only 
trivial ones as for q = 4. For example for q = 6, the terms occuring in the coefficients 
giving k m are even more frequent at order e 4 because of the existing special combinations 
kj + R^/^kj + R 2 . n /^kj = 0. However, in all cases we can write, for k = ki + k' G o"2,i 

Q 

(A e [/)( kl+k ') = ^ 7j (k,e)C/ (k ^ +k ' ) . (44) 
3=1 

Remark 11. The argument e in 7j(k, e) only refers to the perturbation of P 2 CqP 2 in A e 
(see $E3)), and not on the fact that P 2 also depends on e via the radii of the little discs 
composing the set a 2 which are <3(e 1//4 ). 

Remark 12. Due to the form of orders e 2 and e 4 in A € and cf \ and because of (43$ , we 
notice that the dependency in k of the coefficients 7j(k, e) only occurs at orders e° and e . 
Indeed, the dependency in k comes from operators Co at order and {PqCqPq)~ 1 in the 
term P 2 {a(P Q C Q Po)- 1 P Q {aU)} at order e 4 . 

The property that A e U is invariant under the rotation Re. and the identity 

9 

R-nk + kj - ki = R^_ (k + kf +1 - k 2 ) 

lead to 

2 i 

(A t U)M = J^^-k^L^ 1 ^ 1 - 1 ^, k = tf|k G (72,2- 

3=1 

Choosing o~ 2 ^ 2 3 k = k! + k 2 , \k!\ < 8\, we then have 

2g 

(A e [/) {k ' +k2) = y7,(«zik' + ki,f)P (k ' +kj+l) . (45) 
3=1 
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In the same way, after identifying j + 2q with j, we obtain for r = 1, 2q 



= Z)7j+i-r (^(i-o^ + ki.e U7 (k ' 
.7=1 V 9 7 



The important result is that, for a fixed k = ki + k , k G Si, the subspace i?2,k' is invariant 
under the operator A e then denoted . Hence the 2q x 2q matrices of A), are uncoupled 
for different k' G Si. We notice that if jj were independent of k, the lines of the matrix of 

Ae would be deduced each from the previous one by a simple right shift. 

The next useful property of A e is its self-adjointness in E 2 with the Hilbert structure 
of Hq. This property is immediate from the definition (|29p with the scalar product of 



(2) 

the space H s for s = 0. It should be noticed that the full linear operator C e ^ acting in 
E 2 is not selfadjoint in general. Now, isolating the coordinates [/( k ' +kj ), j = l,...,2q for 

( kM 

k' G Si, we still have, for any fixed k , a 2q x 2q self-adjoint matrix A^ due to the previous 
self-adjointness of the operator A e in T~Lq. It results that we have 



lj+l-r \ R «(i-r) k' + ki, e^j = 7 r +l-j \ R *(i-i) k ' + k l, e 



(46) 



We sum up these results in the following 



(S) 

Lemma 13. The subspace E 2 ofH s consisting of functions invariant under rotations by 
n/q may be decomposed into the following Hilbert sum 

E { 2 S) = e E 2jk/ 



k'es 



where we identify the wave vector k with i?ik i.e., ki + k' G a 2) i with k, + R ^u-ij k' G a 2 j 
(see j. 

The 2q— dimensional subspace E 2 y is invariant under the operator A e . Defining the 
coefficients jj(k,e) by 

(Ae[/) (k' +kl ) = ^ 7 . (k)e)C/ (k'+k,) ) k £ a2i 
i=i 

the 2q x 2q matrix of the restriction Ag of A e to i?2.k' *s symmetric and satisfies 
/or any k' G Si. 
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Let us define A Q which is a diagonal matrix with 

7i(k, 0) = (|k| 2 - l) 2 , 72 (k, 0) = 0, .... 72g (k, 0) = 0. 
For k = k' + kj G &2,j, we define 

/3 J (k / ) = (|k' + k J f-l) 2 , j = l,...,2q. 

Hence for j = 1, 2q 



(47) 



7i fiW-j)k' + ki,0j 



Rn(i-j~) k + ki 



1 =/3 J (k') 



(2kf • k' + |k'| 2 ) 2 



fk'") 

and Aq reads 



/ /?i(k') 
/? 2 (k') 



A, 



(k') 


















/3;(k') 











\ 





^-i(k') 




02,00 I 



Then, according to the definition of A^ , we can write 

AC k, )=Af ) +e !, A 1 . 

According to (|46p . in the case when the coefficients 7j(k, e) are independent of k, (which 
corresponds here to the order e 2 ), this leads to a first line for the 2q x 2q matrix, of the 
form 

7u 72> -7 gI lq+i, 7 9 > 7,-1. -73= 72 (48) 

and next lines are deduced by a right shift, making a symmetric matrix. For example in 
the case q = 4, we obtain for Ai a matrix of the form (easily generalizable for any q) 



(11 


72 


73 


74 


75 


74 


73 


72 \ 


72 


7i 


72 


73 


74 


7s 


74 


73 


73 


72 


7i 


72 


73 


74 


75 


74 


74 


73 


72 


7i 


72 


73 


74 


75 


75 


74 


73 


72 


7i 


72 


73 


74 


74 


75 


74 


73 


72 


7i 


72 


73 


73 


74 


75 


74 


73 


72 


7i 


72 


V 72 


73 


74 


75 


74 


73 


72 


7i / 



(49) 



where 7 1; ...7 ?+ i are independent of k'. 
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2.6.2 Computation of coefficients 7- in Ai 

Let us compute the operator 

E ( 2 S) 3 U^e 2 P 2 (aU), 
a = 3t*g — 3(2g— 1). 

where no is given by (HU). We have for U £ (see ([39]) ) 

P 2 ,i(4U) = 2^ 2 ,i(x) + [/ 2 , g+1 (x)e 2ikl - x + 2U 2 ^)e l ^~ k ^ x , 

j=2,..,q,q+2,...,2q 

where we denote by P2 1 the orthogonal projection corresponding to the part 02,1 of the 
spectrum. Since we have 



C/ 2ij (x)e i ( kl - k ^- x = Y U^e 1 ^ 1 -^* = Y u( 

kgo- 2 ,j ketr 



k+kj— ki) ik-x 



2.1 



we obtain 



P 2A (aU) = Y 3eik ' X { U(k) + ^ (k ~ 2kl) + Y 2C/( k+k ^ kl ) \ . (50) 

keCT 2 ,l I j=2,..,q,q+2,...,2q 



As expected, it appears that the linear operator 



U^P 2 (e 2 a)U, U £E ( 2 S) 



(51) 



leaves invariant the subspaces £?2k' and in this subspace it takes the form of a matrix with 
4 identical blocks for the set of 2q coordinates f7( 1,k ), and coefficients jj are independent 
of k' and have the form (1481) with: 



7l=7l+g = 3, 72 =73 = "-7,3 



6. 



(U.') . . . (kM 

Hence, A\ is independent of k' and we write A^ = A±. 

For example, in the case q = 4, we have the following corresponding matrix for 



A, 



-> (P 2 (a- 


)^)U, U e E 2t 


k' 


/I 2 


2 2 12 2 


2 \ 


2 1 


2 2 2 1 2 


2 


2 2 


1 2 2 2 1 


2 


2 2 


2 12 2 2 


1 


1 2 


2 2 12 2 


2 


2 1 


2 2 2 1 2 


2 


2 2 


1 2 2 2 1 


2 


V 2 2 


2 12 2 2 


1 / 



for each fixed k' 6 Si. 
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2.7 Eigenvalues of A e 

From (|29p we have 

AC k, )=A^ + e 2 A 1 . 

We show below the following 

Lemma 14. For any given q > 4, and k € o~2 i/ie eigenvalues of ta&e i/ie /orm 

^ = (2k J -k' + |k , | 2 ) 2 + 3e 2 + 0(e 4 ), 

= [|k' + k/-l] 2 + 3e 2 + 0(e 4 ), j = l,..,2q. 

Proof. The eigenvalues \i € M of Af k ^ satisfy for a certain £ £ M 2? 

{A( k ' ) +6 2 A 1 }C = K- (52) 

Since we deal with selfadjoint operators, any eigenvalue takes the form (see [Kat95| in the 
2g-dimensional subspace E^y.) 

= ^i,o( k ') + eVj,i( k ') + °( e4 )' = "•' 2 ^ 



with 



^• (k') = (2k,-k / + |k'| 2 ) 2 = /3 J (k') 



by definition (j47j) . Eigenvectors take the form 

C^C^o + ^i + Ote 4 ),^!,-^, 

with 

Cj,o = (°) ■•) °' 1) 0.., 0)* , 1 taking the jth place. 
A simple identification at order e 2 leads to 

(AP - Ho )C hl + (A x - Hl )C h0 = 0, (53) 
Taking the scalar product of (|53p with ^ - )0 gives, taking into account the form of Ai, 

_ (AiCj.o.Cj.o) _ 

which is independent of k', and which gives the result of Lemma [T4"l | 

2.8 Inverse of C e in % s 

We already have the following estimate in %q : 

Lemma 15. For any given q > 4, and for e small enough, the linear operator A t is invertible 
in %q with 
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The proof of Lemma [T5l follows directly from Lemma [HI for e < cq, since k' is bounded, 
and all eigenvalues for k' G Si are positive and larger than 2e 2 . 
For extending the estimate to T~L S , we need next property 

Lemma 16. For any K > 0, |x — y\ < K, x and y > 0, and any p > there exists 
d(p, K) > such that 

|(1 + xf - (1 + y) p \ < d(p, K)(l + xf-\ 

The proof of this Lemma is in Appendix A. 
Then we prove the following 

Lemma 17. For any given q > 4, and for e small enough, the linear operator A e is invertible 
in % s for s > 0, with c s > such that 

IIA^II < — 

\\ Il e lis ^ g2 - 

Proof. Let us assume that / G T-L s , and define / s G %o by its Fourier coefficients 

/« = (i + *2)' /2 / (k) - 

We then have ||/J = Then A e U = f leads to {A e U) s G H , ||(A e C/) s || = ||/|| s . 

By definition 

(KU)^ = (1 - |k| 2 ) 2 [/« + e^i^a^C/W, 

where k G 02- Now 

(A e C/)i k ) - [KUs)^ = e 2 $W(l + N 2 ) s ' 2 - (1 + N?)'' 2 ]aP'- r )U®, 
and since |Ak — iVj| < 2 from the form of a, we have from Lemma [TBI 

l(i + a^) s / 2 - (i + N?y' 2 \ < d( s /2, 2)(i + N^y/ 2 - 1 . 

Now define J7 for any U G % s by 

uM = \uM\. 

Then ||[7|| s = \\U\\ B and since for < s < 2, (1 + A 2 ) 5 / 2 " 1 < 1, 

\[(A e U)s ~ (.A € Us)} (k) \ < e 2 d{s/2,2))(aU)W 
(where a differs from the one defined at Lemma [9|), hence for < s < 2 
\\(A e U) s - (A e U s )\\o < d(s/2,2)e 2 \\aU\\ < d s e 2 \ \U\ \ . 

Hence we obtain 

||A e C/ s || < ||(A £ C/) s ||o + d s /2||/|| = ||/|U + d fl /2||/|| , 
and finally, for < s < 2 

\\U\\ s = \\U 3 \\ <^\\A e U s \\ < Jll/IU- 
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Let us prove by induction on s > that ||A e 1 || s < c s e 2 . This holds for < s < 2. Assume 
that it holds for s — 2, then 

|[(A e I7) a - (A e U s )]M\ < e 2 d{s/2,2))(l + Nl) s l 2 - l (aU)M. 

Hence, we have 

\\{KU) S - (A e t7 s )|| < d{s/2,2)e 2 \\aU\\ s . 2 < d s e 2 \ \U\ | s _ 2 . 

We assumed that 1 1 1 1 1 s 2 < ^7^, hence we obtain 

||A e t/ s ||o < ||(A e [/) s || + 4c s - 2 ||/|| s -2 = H/IU + d s c s _ 2 ||/|| s _ 2 , 

hence 



U s \\o < -^\\A e U s \\o < - 2 
This ends the proof of Lemma H71 | 
Then we finally have 

Lemma 18. For any q > 4, and s > 0,there exists eo > 0, such that for < e < eo the 
linear operator C t has a bounded inverse in fi s , with 

ll^e lis ^ £ 2 ' 

where c(s) is a positive constant only depending on s. 
Proof. From Lemma [T7] we have 

(£®)- 1 = 0- + e*K 1 Rt)- 1 K\ 

and ||e 4 A~ 1 i2 e ||<, < c s e 2 ||i? e || s < c' s e 2 . For e small enough we then have 

wm-'w < ^. 

Then, from Lemma [9] we deduce immediately that there exists a constant c(s) such that 

<>), 



which proves the Lemma. 



3 Existence of the solution 

Below we prove our main result 

Theorem 19. For any q > 4 and for any s > q/2, there exists Xq > such that for 
< A < Ao, there exists a quasipattern solution of the Swift- Hohenberg steady equation 
07]) in % s , bifurcating from and invariant under rotations of angle ir/q. Its asymptotic 
expansion at the origin is given by the formal expansion computed in JIRlOj/ . 
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Proof. We want to solve (|2ip with respect to W in H s ,s > q/2. Taking into account of 
Lemma [T8| this equation takes the following form for e ^ : 

W = -e 3 C~ l [f e + 3eU e W 2 + e 5 W 3 } , (54) 

which we write as 

W = Q(e,W) 

where Q is well defined in ]0, eo[x7i s , depending smoothly on its arguments for < e < eo 
and W E H s . In fact for e ^ fixed, Q is analytic in W, and observing that e^Cj 1 = 0(e), 
we see that Q is continuous in e on ]0, eo[. 

The map Q is Lipschitz in W in a fixed ball of T-L s , with a small Lipschitz constant for 
e e]0, eo[. Indeed, we have 

\\G(e, W) - G(e, W')\\ s < W^C-^eU^W 2 - W 2 ) + e 5 (W 3 - W ,3 )]\\ s , 

e 3 ||£~ || s < c s e and \\U e \\ s < c' s e. Moreover, we have 

g( e ,0) = -e 3 C 1 (A) = 0(6). 

Than according Dieudonne's version of the implicit function theorem [D60] [(10.1.1)], 
there exists a unique mapping W(e) into a ball in % s of size O(e), such that 

W(e) = G(e,W(e)) 

for all e E]0, eo[ and W is continuous there. Finally we have a solution U = U e + e 4 VF(e) of 
(fTTD of the form ((20]). | 



A Proof of Lemma [6] 

We follow and modify the argument of Iooss-Rucklidge [IR10] [appendix C] used to prove 
that % s is an algebra jlRlO] [lemma 4.2]. Let 



ker 



e ik.x v 



5> (k) H 

ker 



,ik.x 



be elements of Tl s n Tl s i. We have 



2- 2s+1 ||™||2<^ 



K 



(k)„(k') 



k+k'=K 



(i+jvg)-+x; 



K 



E 

k+k'=K 



n (k) v (k') 



(l+iv^r- 



Si 



5, 



Moreover, we have 



K 



E 



k+k'=K 

W k <3iV k , 



K 



E 



k+k'=K 

iV k >3JV k , 



(i + N 2 y. 



s> 
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Using the fact that 1 < 



/ 16(l+JVg,) 
I (1+^k) 



and Cauchy-Schwarz inequality, we obtain (s' > q/2) 



K 



16 



< C||u|| 2 |M| 2 ,. 



To obtain a similar bound for S'{, we use the Iooss-Rucklidge dyadic decomposition of 
S" : A p u := ^ 2 p<Af k <2P+ 1 u^e^ x , A_iu := and SkU := Ylp=-i ^p u - Then, u = 
Y, P >-i A p u G Hs if and only if Y P >-i 2 2ps ||A p -u||q < +00. According to the computation 
of pin] [p. 387], we have 



( 



E 

K 



E 



u (kyk') 



■ k+k'=K 

\Ar k >3Af k , 



JK.x 



/ 



< C 2 2J ' 
i=-i 



j'+i 
. p=j-i 



(55) 



Since s' > q/2, by Cauchy-Schwarz, we have Yl M k ^\ < C IMU'- So, following the computa- 
tions [IR10] [p. 388], we obtain [|/Sp_i«A p it||o < C||Apit||g[|u||^/. Prom this and following the 
same computation as in [TRIO] [p. 388], we obtain 



S'(<C"\\v\\ 2 s , 2 2ps ||A p n||2 < Ci| 
p=-l 



n2 11 n2 
V\\ .1 \\u\\ - 



To get an estimate for 52, we just need to interchange the role of u and v and the result is 
proved. 



B Proof of Lemma [16 

We assume x, y > 0, and p > and 

\x — y\ < K. 

Then, we prove that 

\{l + xf -{l + yf\<d{p){l + x) p -\ 

with 

d(p) = pK(l + K) for p > 1, 

= pK(l + K) 1 " p forp< 1. 

Proof. For some t between x and y, we have 

(1 + x) p - (1 + y) p = p(x -y)(l + tf-\ 

If p > 1 we use (1 + t)/(l + x) = 1 + (t - x)/(l + x) < 1 + K if x < t < y and < 1 if 
y < t < x. This proves the lemma if p > 1. 

Similarly if p < 1 we use (l+x)/(l+i) < 1 if x < f < y and = l + (x-t)/(l + t) < l + K 
if y < t < x and so (1 + t)^ 1 < const. (1 + x) p_1 . | 
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